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Small radii of neutron stars as an indication of novel in-medium effects 


Wei-Zhou Jiang/’^’* Bao-An Li,^A and F. J. Fattoyev^’^A 
^Department of Physics, Southeast University, Nanjing 211189, China 
^Department of Physics and Astronomy, Texas A&M University-Commerce, Commerce, TX 75j29, USA 
^Center for Exploration of Energy and Matter, Indiana University, Bloomington, IN j7f08, USA 

(Dated: September 14, 2015) 

At present, neutron star radii from both observations and model predictions remain very uncertain. 
Whereas different models can predict a wide range of neutron star radii, it is not possible for most 
models to predict radii that are smaller than about 10 km, thus if such small radii are established 
in the future they will be very difficult to reconcile with model estimates. By invoking a new term 
in the equation of state that enhances the energy density, but leaves the pressure unchanged we 
simulate the current uncertainty in the neutron star radii. This new term can be possibly due 
to the exchange of the weakly interacting light U-boson with appropriate in-medium parameters, 
which does not compromise the success of the conventional nuclear models. The validity of this new 
scheme will be tested eventually by more precise measurements of neutron star radii. 

PACS numbers: 26.60.-c, 14.70.Pw, 97.60. Jd 


Neutron star (NS) is a unique place to test fundamen¬ 
tal forces at the extremes of matter density, gravity and 
magnetic fields. Unfortunately, uncertainties in both the 
Equation of State (EOS) of super-dense nuclear matter 
and the strong-held gravity strongly interplay with each 
other in determining observational properties of neutron 
stars, for the latest review, see, e.g., [1]. For instance, 
in a simple version of modihed gravity where the non- 
Newtonian gravity exists, neutron stars could have very 
different structures compared to predictions using Ein¬ 
stein’s General Relativity (GR) theory of gravity [2-4]. 
The radius of a neutron star is one of the most important 
observables sensitive to the underlying nuclear EOS and 
gravity theories used. Gurrently, within GR the radius 
of a canonical NS has been predicted to be roughly from 
11 to 15 km [5-9] depending on the EOS used. Provided 
the third family of compact stars known as strange stars 
exist, their radii could be as small as 7 or 8 km [10- 
12], although these models normally predict star masses 
much smaller than the masses of observed massive neu¬ 
tron stars. Thus, the measurement of NS radii plays 
a very important role in resolving several issues in fun¬ 
damental physics. Unfortunately, the extraction of NS 
radii from observations still suffers from large systematic 
uncertainties [13] involved in the distance measurements 
and theoretical analyses of the light spectrum [6, 14-16]. 
Gonsequently, a wide range of the radius with the mass 
around 1.4 Mq has been reported [16-24]. In particu¬ 
lar, using the thermal spectra from quiescent low-mass 
X-ray binaries (qLMXBs) Guillot and collaborators ex¬ 
tracted NS radii of i?NS = 9.4 ± 1.2 km [24]. Another 
recent comprehensive study of spectroscopic radius mea¬ 
surements suggest that for a 1.5 Mq NS the extracted 
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radii are 10.® km [25]. It is important to note that 
at the moment no consensus has been reached yet on the 
extracted NS radii. For instance, Bogdanov found a 3-cr 
lower limit of 11.1 km on the radius of the PSR J0437- 
4715 [26], and Poutanen et al. got a lower limit of 13 
km for 4U 1608-52 [27]. Whether the radii of canoni¬ 
cal neutron stars can be as low as 10 km have been an 
interesting topic of hot debate during the last few years. 

Whereas the situation has been signihcantly improved 
over the last few years, the systematic errors have been 
hindering severely the accurate determination of the NS 
radii from astrophysical observations. However, if the ex¬ 
istence of NSs with small radii is firmly established, they 
would pose a severe challenge to the current models of 
the nuclear EOSs. While it is not very difficult to satisfy 
the maximum mass constraint, to our best knowledge, no 
nuclear models available with best-fit parameters to date 
can reproduce the small NS radius constraint. There are 
only a few microscopic models or approaches that— ei¬ 
ther disregard some of the nuclear physics constraints [5] 
or adjust the high-density part of the EOS using various 
polytropes to match the whole density profile obtained 
from observation [28]— can account for both constraints. 
From the nuclear physics standpoint, the small NS ra¬ 
dius requires certain softness of the EOS of the isospin- 
asymmetric nuclear matter. As one of the basic blocks 
of the EOS of asymmetric matter, nuclear symmetry en¬ 
ergy around 1-2 times the saturation density of nuclear 
matter plays a dominating role in determining the radii 
of neutron stars [6]. A softening of the symmetry en¬ 
ergy can lead to an appreciable decrease of the NS ra¬ 
dius [7]. However, with the maximum NS mass held ap¬ 
proximately at a constant, most non-relativistic and rel¬ 
ativistic models that are facilitated with soft symmetry 
energies can only bring the radii of canonical NSs down 
to about 12-13 km [8, 9, 17, 21]. Moreover, it becomes 
very difficult to further reduce the NS radius by further 
softening the symmetry energy. In particular, one would 
then encounter the stability problem in the NS matter 
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EOS when the symmetry energy becomes too soft [3]. 
To further reduce the NS radius, one could imagine to 
reduce the pressure of the isospin-symmetric part of the 
EOS in the intermediate density region. But the space 
in so doing is actually limited by the saturation prop¬ 
erties of nuclear matter and the constraint on the EOS 
of dense nuclear matter extracted from studying nuclear 
collective flow in high energy heavy-ion reactions [29]. 
Furthermore, as we can see from the empirical relation 
Rp~^/'^{Pb) ~ C{pb) between the NS radii R and the 
pressure p{pb) with C{pb) being a constant at a given 
baryon density pB [7], such a reduction is also rather in¬ 
efficient, since the isospin-symmetric EOS contribution 
to the total NS matter pressure p{pb) is relatively small 
in the relevant density region, where this empirical re¬ 
lation holds. Moreover, even if the signihcant reduction 
were allowed for the total pressure, the limited decrease 
of the NS radius would be at the cost of a large reduction 
of the NS maximum mass, because the significantly re¬ 
duced pressure needs the corresponding reduction of the 
NS mass to balance the gravity therein, and also because 
the NS mass is the total energy integrated in a nutshell 
with the reduced radius. A significant reduction in the 
maximum mass is certainly disfavored by the recently 
discovered massive neutron stars of 2Mq that require a 
stiff EOS [30, 31]. 

Facing currently with this severe theoretical issue, we 
explore in this work a new possibility to soften the EOS: 
adding a new term in the EOS that enhances the en¬ 
ergy density while keeping the pressure unchanged. In so 
doing, the enhancement of the energy density requires 
some shrinkage in NS radius without reducing signifi¬ 
cantly the NS maximum mass. Since the pressure is 
p = p\d((.j pb)I dPB■, to keep the pressure unchanged we 
modify the energy density e by 

e = Co + ClPb-, (1) 

where eo is the base energy density given by any model. 
The second term linear in density is an addition to the 
base EOS with Cl being the amending coefEcient. As 
we shall explain, this modihcation to the EOS of isospin- 
asymmetric nuclear matter can be realized by considering 
the interaction added by a vector boson with appropriate 
in-medium parameters. 

In principle, the modification given in Eq. (1) could be 
added to any nuclear model. In this work, we just demon¬ 
strate the effects using several typical relativistic mean- 
field (RMF) models. The RMF models under considera¬ 
tion include the SLC, SLCd [32], SL3 [9] and NL3040 [33]. 
The SL3 and the NL3040 have similarly stiff EOSs at 
high density and both give large NS maximum masses of 
more than 2.6Mq, while the SLC and the SLCd feature 
the same EOS of symmetric matter within the constraints 
obtained from analyzing the collective flow in relativistic 
heavy-ion collisions [29] . The only difference between the 
SLC and the SLCd is that the latter has a softer sym¬ 
metry energy. In a similar fashion, the NL3040 was also 
built from the original NL3 to feature a softer symmetry 
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FIG. 1: The relation between the matter pressure and the 
energy density with RMF models SLCd and SL3. Note that 
the displacement of energy density with a given Cn is the 
same at a given density but not at a constant pressure. 


energy. The stiffness of the symmetry energy is normally 
measured by its density slope at the saturation density 
of nuclear matter po, L = 3po i9EsymipB)/dpB)p^- The 
value of L for the SLC and SL3 is 92.3 and 97.1 MeV, 
while it is 61.5 and 45.0 MeV for the SLCd and NL3040, 
respectively. For a comparison, it is interesting to note 
that currently the most probable value of L is in the 
range of 40 < L < 70 MeV according to recent analy¬ 
ses of various terrestrial experiments and astrophysical 
observations, see, e.g., Refs. [34-40] and Ref. [41] for a 
comprehensive review. Thus, the SLC and SL3 are obvi¬ 
ously too stiff while the SLCd and NL3040 are consistent 
with the existing constraints in terms of their L values. 
Nevertheless, they are all appropriate for the purposes of 
this study. 

Shown in Fig. 1 are two examples of the EOS, i.e., 
pressure versus energy density, with the SLCd and SL3 
parameter sets. It is seen that at a constant pressure, 
the amending term can soften the EOS considerably, 
i.e., reducing the slope of the pressure with respect to 
the energy density, especially with the SLCd. However, 
the relative effect of the amending term goes down with 
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the increasing density because it is just linear in density 
while the EOS of usual nuclear models evolves generally 
with the density squared. We emphasize that the EOS 
softening scheme considered here is quite different from 
the usual mechanisms mentioned in the introduction. In 
particular, typical phase transitions to matter with new 
degrees of freedom normally reduce the maximum mass 
dramatically, but often keep the NS radius more or less 
the same because the phase transitions usually occur in 
the small inner core of NSs. Of course, exceptions may 
exist when the new degrees of freedom, such as the A 
resonances, can emerge at a relatively low density [42]. 
It is interesting to note that by using Lindblom’s inver¬ 
sion algorithm [43] Chen and Piekarewicz were recently 
able to obtain a softened EOS from the given small NS 
radii [44]. 

We now examine effects of the amending term on the 
radii of neutron stars. As in Ref. [20], here we consider 
the simplest model of neutron stars consisting of just 
neutrons, protons and electrons. Shown in Fig. 2 are the 
mass-radius (MR) trajectories of neutron stars with the 
amending term within various RMF models. The amend¬ 
ing coefficient is exemplified as Cl = 100, 200 and 300 
MeV, and the results with original models are displayed 
with Cl = 0. Comparing results of the original models, 
we see that the softening of the symmetry energy may re¬ 
duce the NS radius by as large as 1.5 km for a canonical 
NS when L is reduced from 97.1 to 45 MeV. The space 
for further reducing the slope parameter L is small, and 
in fact, the further reduction in L has a very limited ef¬ 
fect in decreasing the NS radius. Moreover, it is seen 
that even with the significant softening of the symmetry 
energy within the original models the NS radius is still 
far above 9.4±1.2 km extracted by Guillot et al. [24]. If 
such small radii are established, it is then interesting to 
see that the amending term can indeed further reduce the 
NS radius. Obviously, the role of the amending term is 
similar in all models: the larger the amending coefficient 
is, the more is the reduction of the NS radius. With the 
same amending coefficient in different models, the shifted 
magnitude is also similar. Typically, the amending coef¬ 
ficient, varying up to 300 MeV, can cause a reduction of 
about 3 km in the NS radius. With the amending coeffi¬ 
cient of Cl = 300 MeV, we see that the MR trajectories 
with the SLC and SLCd fall into the regime extracted 
by Guillot et al. [20]. We see from Fig. 2 that the NS 
maximum mass is still not reached in the SLC and SLCd 
models before reaching the causality boundary, because 
the allowed nucleon density has a maximum in the con¬ 
struction of such models to meet the chiral limit [9] . The 
removal of such a limiting density may bring the NS max¬ 
imum mass closer to the causality limit. With the larger 
amending coefficient, the MR trajectories for the SL3 
and NL3040 can be very close to the upper margin of the 
extracted regime. Notably, we see that the MR trajec¬ 
tories are not clearly away from the causality constraint, 
though the amending coefficient causes the decrease of 
the NS maximum mass. Here, the moderate reduction of 


the NS maximum mass is just because of the softening 
of the EOS, namely no excess of pressure can resist the 
additional gravity arising from the increase of the energy 
density. 
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FIG. 2: (Color online) The mass-radius trajectories of neu¬ 
tron stars with RMF models: SLC, SLCd, SL3 and NL3040. 
In each panel, the amending coefficient Cl is taken to be the 
values 0, 100, 200, 300 MeV, respectively. The slope param¬ 
eter of the symmetry energy is also labeled in each panel. 

In the RMF framework, the amending term in Eq. (1) 
can be understood as a specific in-medium effect. Similar 
to the analysis in Refs. [45-47], the amending term, in¬ 
corporated into the vector potential, leads to the density- 
dependent coupling constant of the vector (w) meson 

gliPB) = {gujU} + 2CL)ml/pB. ( 2 ) 

This relation indicates that the larger Cl is responsible 
for the stronger density dependence of the coupling con¬ 
stant, and with the increase of density, the in-medium 
effect decreases with the growing uj. Fig. 3 shows the 
density-dependent coupling constant for the models con¬ 
sidered. For a comparison, the gojips) obtained from the 
Dirac-Brueckner (DB) potential of Bonn A is also shown 
in the figure. It is seen that the density dependence, simi¬ 
lar to the one from the DB potential, is needed to produce 
a significant reduction of the NS radius. We can infer, in¬ 
deed, that the density dependence here is stronger than 
that from the DB potential, because the latter owns a 
partial cancelation of the density dependencies between 
the scalar and vector potentials [45-47]. In addition, we 
see that the in-medium effect of all cases tend to vanish 
at high densities. This means that the decrease of the NS 
radius is dominated by the modification of the amending 
term to the low-density EOS. For the dropping of the NS 
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maximum mass, it can then hopefully be cured by mod¬ 
ifying the high-density component of the EOS of nuclear 
matter. 
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FIG. 3: (Color online) The density-dependent vector coupling 
constant incorporated from the amending term in RMF mod¬ 
els: SLC/SLCd and NL3040. The density-dependent coupling 
constant obtained from the DB potential of Bonn A is also 
given for comparison. 


The DB potential is obtained by solving the two-body 
correlations. One may expect that the many-body rather 
than just two-body correlations may generate more in¬ 
medium effects [5]. It should, however, be pointed out 
that the arbitrarily strong density dependence usually 
lacks observable grounds. In the present RMF models, 
one should evade the violation of the low-energy con¬ 
straints in finite nuclei from the in-medium effect induced 
by the amending term. One possible solution is to resort 
to the exchange of the very light and weakly interacting 
boson that is undetected to date. A favorite candidate is 
the U-boson, which is actually a particle beyond the stan¬ 
dard model [48, 49]. The light U-boson, first proposed by 
Fayet [48], might be regarded as the mediator of the pu¬ 
tative fifth force [49-51]. Recently, the possibility of the 
MeV dark matter with the light U-boson mediation was 
considered to account for the bright 511 keV 7 emission 
of positron annihilations from the galactic bulge [52-57] , 


despite that there are a number of conventional astro- 
physical sources for positron annihilations, see Ref. [58] 
and references therein. In the past, the signihcant effects 
of the U-boson in neutron stars were predicted [2-4], al¬ 
beit the boson should couple to baryons very weakly. The 
modihcation of the EOS due to the U-boson in the mean- 
held approximation simply reads 


1 



(3) 


with gu/rUu being the ratio of the U-boson coupling con¬ 
stant and its mass. 

By assuming a density-dependent U-boson mass of 
we realize the linear density depen¬ 
dence of the amending energy density. If the U-boson 
has a very light bare mass of m„o in free space, its in¬ 
medium mass can be given as 

ml = glpB/2CL + mlo- (4) 


In general, the in-medium effect of the boson parame¬ 
ters is associated with the contribution of the intermedi¬ 
ate states. In the RMF theory, the in-medium effect is 
usually attributed to the nonlinearity of the meson self¬ 
interactions [59]. The in-medium effect for the boson 
may also be realized by carefully choosing the nonlinear 
self-interacting terms. We leave this problem for a future 
study. 

Regarding the possible violation of the low-energy nu¬ 
clear constraints for hnite nuclei, one can avoid it by lim¬ 
iting the weak interaction strength of the U-boson. For 
instance, if gu is 0 . 01 , the interaction strength is about 
two orders of magnitude weaker than the electromagnetic 
interaction, and such weak interaction is not able to af¬ 
fect properties of finite nuclei. In this sense, the present 
scheme to invoke the U-boson does not compromise the 
success of optional nuclear models. The interesting ques¬ 
tion is what behavior of the U-boson mass will allow 
such a weak interaction strength. Shown in Fig. 4 is 
the density-dependent mass of the U-boson as a function 
of density, for = 0.01. Here, as an example, the small 
bare boson mass is taken to be as 0.1 MeV, which is just a 
free parameter. We see that for the given coefficients Cl, 
the in-medium mass is just within a few MeV. With the 
increase of the Cl, the in-medium mass of the U-boson 
decreases. For Cl = 300 MeV, the U-boson mass is close 
to 1 MeV, which is consistent with that considered in 
Ref. [4]. 

We should say that the role of the vector boson in soft¬ 
ening the EOS eventually is rather pragmatic, albeit pu¬ 
tatively accounting for small NS radii. Usually, the vector 
boson is a source of the pressure, while here it gives zero 
contribution. From the point of view of nuclear many- 
body approaches, the repulsion of the vector meson can 
be softened in the intermediate density region due to the 
correlation effect from the intermediate-state contribu¬ 
tions, and such a softening can be simulated in the RMF 
theory by invoking the nonlinear self-interacting terms. 
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At high densities where the intermediate-state contribu¬ 
tion is small due to the Pauli blocking the reduction of 
the in-medium effect of the vector boson will then recover 
the repulsion which may enable a stiffer EOS at high den¬ 
sities. In this sense, if we use a more complicated density- 
dependent term in Eq.(l) in the high density region, the 
corresponding U-boson can then stiffen the high-density 
EOS and the U-boson mass would be also much closer 
to a constant at high densities. While we use a simple 
amending term in Eq.(l) to demonstrate the reduction 
of the NS radius, it should not represent that the corre¬ 
sponding vector U-boson would simply soften the EOS. 



FIG. 4: (Color online) The in-medium mass of the U-boson 
with various Cl as a function of density. Here, Qu = 0.01 and 
muo = O.lMeV. 

The NS radius can attain larger reduction by further 
increasing the amending coefficient. In our analysis here 
we have invoked the special amending term that does 
not affect the pressure. Generally, other forms could also 
be optional, as long as they result in the reduction of 
the NS radius. A favorable form should reduce the low- 
density pressure and satisfy all constraints on the EOS 
of nuclear matter at saturation density. To meet these 
demands, one can consider the light scalar boson. Note 
that the scalar boson, like the scalar meson, can also 


modify the mass of baryons. We have checked that the 
inclusion of such a light scalar boson can again reduce 
the NS radius by 0.5-1.5 km, depending on the model 
used. To keep the pressure positive definite, the coupling 
constant of the scalar boson should be much weaker (e.g., 
gs = 0.005) than that of the vector boson, and the NS 
maximum mass is little changed by the scalar boson. 

In summary, we have demonstrated that the large un¬ 
certainty of the NS radius can be simulated by introduc¬ 
ing the amending term that enhances the energy density, 
but leaves the pressure unchanged. The distinct effect 
of the amending term on the reduction of the NS radius 
is rather universal to any nuclear EOS. The new term is 
linear in density, and the incorporation into the relativis¬ 
tic mean-field potential produces the strong in-medium 
effect on the vector coupling constant at low densities. 
The strong in-medium effect can be responsible for the 
significant reduction of neutron star radii, and may pro¬ 
vide an explanation for the small NS radius if it is firmly 
established. We interpret this novel in-medium effect by 
the possible exchange of a putative weakly interacting U- 
boson. By requiring that the low-energy constraints in 
finite nuclei must not be violated, we find interestingly 
that the in-medium U-boson mass is just below a few 
MeV with the interaction strength being two orders of 
magnitude weaker than the electromagnetic interaction. 
Finally, the validity of the scheme presented here will be 
tested eventually by the more precise measurements of 
the NS radii. 


Acknowledgement 

We would like to thank W.G. Newton for helpful dis¬ 
cussions. The work was supported in part by the Na¬ 
tional Natural Science Foundation of Ghina under Grant 
No. 11275048 and 11320101004, the China Jiangsu 
Provincial Natural Science Foundation under Grant No. 
BK20131286, the US National Science Foundation un¬ 
der grants No. PHY-1068022 and PHY-1205019, the 
U.S. Department of Energy’s Office of Science under 
Award Number DE-SC0013702, the CUSTIPEN (China- 
U.S. Theory Institute for Physics with Exotic Nuclei) 
under DOE grant number DE-FG02-13ER42025, and 
DOE grants DE-FG02-87ER40365 (Indiana University) 
and DE-SC0008808 (NUCLEI SciDAC Collaboration). 


[1] X.T. He, F. J. Fattoyev, B.A. Li and W. G. Newton, 
Phys. Rev. C 91, 015810 (2015). 

[2] M. I. Krivoruchenko, F. Simkovic, and A. Faessler, Phys. 
Rev. D 79, 125023 (2009). 

[3] D. H. Wen, B. A. Li, and L. W. Ghen, Phys. Rev. Lett. 
103, 211102 (2009). 

[4] D. R. Zhang, P. L. Yin, W. Wang, Q. C. Wang, and W. 
Z. Jiang, Phys. Rev. G 83, 035801 (2011). 


[5] R. B.Wiringa, V. Fiks and A. Fabrocini, Phys. Rev. C 
38, 1010 (1988). 

[6] J. M. Lattimer and M. Prakash, Astrophys. J, 550, 426 

( 2001 ). 

[7] J. M. Lattimer and M. Prakash, Phys. Rep. 442, 109 
(2007). 

[8] B. A. Li and A. W. Steiner, Phys. Lett. B642, 436 (2006). 

[9] W. Z. Jiang, B. A. Li and L. W. Ghen, Phys. Lett. B653, 







6 


184 (2007). 

[10] B. Kampfer, Phys. Lett. BlOl, 366 (1981). 

[11] N. K. Glendenning and C. Kettner, Astron. Astrophys. 
353, L9 (2000). 

[12] J. Schaffner-Bielich, M. Hanauske, H. Stocker, and C. 
Greiner, Phys. Rev. Lett. 89, 171101 (2002). 

[13] M. C. Miller, arXiv/1312.0029. 

[14] P. Haensel, Astron. Astrophy. 380, 186 (2001). 

[15] G. M. Zhang, H. X. Yin, Y. Kojima, H. K. Chang, et.ah, 
Mon. Not. Roy. Astron. Soc. 374, 232 (2007). 

[16] V. Suleimanov, J. Poutanen, M. Revnivtsev, and K. 
Werner, Astrophys. J. 742, 122 (2011). 

[17] A. W. Steiner, J. M. Lattimer, and E. F. Brown, Astro¬ 
phys. J. 722, 33 (2010). 

[18] F. Ozel, G. Baym, and T. Guver, Phys. Rev. D 82, 
101301 (2010). 

[19] T. Giiver and F. Ozel, Astrophys. J. 765, LI (2013). 

[20] S. Guillot, M. Servillat, N. A. Webb, R. E. Rutledge, 
Astrophys. J. 772, 7 (2013). 

[21] J. M. Lattimer and A. W. Steiner, Astrophys. J. 784, 
123 (2014). 

[22] A. W. Steiner, J. M. Lattimer, and E. F. Brown, Astro¬ 
phys. J. Lett. 765, L5 (2013) 

[23] C. O. Heinke, H. N. Cohn, P. M. Lugger, et. ah, Mon. 
Not. Roy. Astron. Soc. 444, 443 (2014). 

[24] S. Guillot and R. E. Rutledge, Astrophys. J. 796, 1, L3 
(2014). 

[25] F. Ozel, D. Psaltis, T. Giiver, G. Baym, C. Heinke, and 
S. Guillot, arXiv: 1505.05155. 

[26] S. Bogdanov, Astrophys. J. 762, 96 (2013). 

[27] J. Poutanen, J. Nattila, J. J. E. Kajava, et. ah, Mon. 
Not. Roy. Astron. Soc. 442, 3777 (2014). 

[28] K. Hebeler, J. M. Lattimer, C. J. Pethick, and A. 
Schwenk, Astrophys. J. 733, 11 (2013). 

[29] P. Danielewicz, R. Lacey, and W. G. Lynch, Science 298, 
1592 (2002). 

[30] P. B. Demorest,T. Pennucci, S. M. Ransom, M. S. E. 
Roberts, and J. W. T. Hessels, Nature 467, 1081 (2010). 

[31] J. Antoniadis, P. C. C. Freire, N. Wex, et al. Science, 
340, 448 (2013). 

[32] W. Z. Jiang, B. A. Li, and L. W. Chen, Phys. Rev. C 76, 
054314 (2007). 

[33] J. Piekarewicz, Phys. Rev. C 69, 041301(R) (2004). 

[34] M. B. Tsang, J. R. Stone, F. Camera, et al. Phys. Rev. 
C 86, 015803 (2012). 


[35] W. G. Newton, M. Gearheart, and B. A. Li, Astrophys. 
J. Suppl. Ser. 204 , 9 (2013). 

[36] W. G. Newton and B. A. Li, Phys. Rev. C 80, 065809 
(2009). 

[37] A. W. Steiner, and S. Gandolfi, Phys. Rev. Lett., 108 , 
081102 (2012). 

[38] D. H. Wen, W. G. Newton, and B. A. Li, Phys. Rev. C 
85 , 025801 (2012). 

[39] J. M. Lattimer, Annu. Rev. Nucl. Part. Sci. 62 , 485 
( 2012 ). 

[40] B.A. Li and X. Han, Phys. Lett. B 727, 276 (2013). 

[41] B.-A. Li, A. Ramos, G. Verde, and 1. Vidana, eds., ’’Top¬ 
ical issue on nuclear symmetry energy”, Eur. Phys. J. A 
50, No. 2, (2014). 

[42] B.-J. Cai, F. J. Fattoyev, B.-A. Li, and W. G. Newton, 
Phys. Rev. C92, 015802 (2015). 

[43] L. Lindblom, Astrophys. J. 398, 569 (1992). 

[44] W.-C. Chen and J. Piekarewicz, arXiv: 1505.07436. 

[45] R. Brockmann and H. Toki, Phys. Rev. Lett. 68, 3408 
(1992). 

[46] R. Fritz, H. Miither and R. Machleidt, Phys. Rev. Lett. 
71 , 46 (1993). 

[47] Z. Y. Ma, H. L. Shi and B. Q. Chen, Phys. Rev. C 50, 
3170 (1994). 

[48] P. Fayet, Phys. Lett. B 95, 285 (1980). 

[49] P. Fayet, Phys. Lett. B 172 , 363 (1986); Nucl. Phys. B 
347 , 743 (1990). 

[50] E. Fischbach and C. L. Talmadge, The Search for Non- 
Newtonian Gravity (Springer-Verlag, Inc., New York, 
1999), ISBN 0-387-98490-9. 

[51] E. G. Adelberger et ah, Annu. Rev. Nucl. Part. Sci. 53, 
77 (2003). 

[52] C. Boehm, D. Hooper, J. Silk, M. Casse, and J. Paul, 
Phys. Rev. Lett. 92, 101301 (2004). 

[53] C. Boehm and P. Fayet, Nucl. Phys. B 683, 219 (2004). 

[54] N. Borodatchenkova, D. Choudhury, and M. Drees, Phys. 
Rev. Lett. 96 , 141802 (2006). 

[55] S. H. Zhu, Phys. Rev. D 75, 115004 (2007). 

[56] P. Fayet, Phys. Rev. D 75, 115017 (2007). 

[57] P. Jean et al., Astron. Astrophys. 407 , L55 (2003); J. 
Knodlseder et ah, Astron. Astrophys. 411 , L457 (2003). 

[58] N. Prantzos, C. Boehm, A.M. Bykov, et.al.. Rev. Mod. 
Phys. 83 , 1001 (2011). 

[59] J. Boguta and A. R. Bodmer, Nucl. Phys. A 292 , 423 
(1977). 


